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Abstract 

We investigate the cosmological and the local tests of the f{R) theory of modified gravity via the 
observations of (1) the cosmic expansion and (2) the cosmic structures and via (3) the solar-system 
experiments. To fit the possible cosmic expansion histories under consideration, for each of them 
we reconstruct f(R), known as "designer /(i?)". We then test the designer f(R) via the cosmic- 
structure constraints on the metric perturbation ratio ^/<1> and the effective gravitational coupling 
G e Q and via the solar-system constraints on the Brans-Dicke theory with the chameleon mechanism. 
We find that among the designer f(R) models specified by the CPL effective equation of state w e s , 
only the model closely mimicking general relativity with a cosmological constant (ACDM) can 
survive all the tests. Accordingly, these tests rule out the frequently studied w e g = — 1 designer 
f(R) models which are distinct in cosmic structures from ACDM. When considering only the 
cosmological tests, we find that the surviving designer f(R) models, although exist for a variety of 
w e s, entail fine-tuning. 
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I. INTRODUCTION 



Our universe is presently at the stage of accelerating expansion, as suggested by the 



cosmological observations 



j, |2j . Assuming a homogeneous and isotropic universe and based 
on general relativity (GR), one cannot explain this late-time cosmic acceleration if the 
universe contains only matter and radiation. A substance with negative pressure adding to 
the right-hand side of the Einstein field equations, known as dark energy [31], can lead to the 
cosmic acceleration |4|. Observations [5H7( suggest that our universe should be nearly flat 
and consist of 73% dark energy, 22% dark matter, 5% ordinary matter, and approximated 
0.008% radiation. It is ironic that so far our knowledge of physics can only explain 5% 
of the universe, with nearly 95% of the universe remained unexplained. This is the most 
challenging issue for the physicists of the 21st century. 

Instead of invoking dark energy, there have been attempts to modify the left-hand side 
of the Einstein ennation^e.g. vian reposing alternative gravity theories, to explain the 



late-time acceleration 



llj (see [12) for a review). These attempts are known as modified 



gravity. Among these theories, in the present paper we will focus on the f{R) theory 
of modified gravity where the Einstein-Hilbert Lagrangian density is added by a general 
function of the Ricci scalar, f{R). 

A. A. Starobinsky was the first to point out that de Sitter space, and therefore the 
exponential expansion of spacetime, can be achieved in the early universe by quantum 
correction, (R + a 2 R 2 ), to the Einstein-Hilbert action [k^. This UV modification leads to 
an early accelerating universe (Starobinsky inflation) and motivates people to construct the 
f(R) models [l5j to explain the late-time acceleration of the universe. Since the f (R) 
theory is equivalent to the Brans-Dicke theory with the Brans-Dicke parameter cubd = 16], 
it seems that this is already ruled out by the constraint wbd > 4 x 10 4 set by observations 
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19J. Nevertheless, it was pointed out recently that in a non- minimally coupled scalar- 



tensor theory the scalar fie 



interaction with matter 
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d may have a scale- or environment-dependent mass via its 



2l| . That is, the higher (lower) the matter density is, the more 



(less) massive the scalar field is. Consequently, in the solar system, which has a much higher 
density than the universe, the scalar field would have a large effective mass (and therefore a 
small Compton wavelength) so that it can evade the solar-system tests. At the mean time, 
the same scalar field would have a much smaller effective mass at cosmological scales, as 
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small as the present Hubble expansion rate (i.e. a long Compton wavelength comparable to 
the size of the universe). Accordingly, it can be a candidate driving the late-time cosmic 



acceleration. The sensitivity to 
name, "chameleon mechanism" 



;he environment of such a scalar field interaction gives the 
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2l|. 



In this paper we investigate various tests of the f(R) theory of modified gravity. As 
an alternative to dark energy for driving the late-time cosmic acceleration, the f(R) theory 
needs to pass the cosmo logical tests, including the constraints about (1) the cosmic expansion 
and (2) the cosmic structure formation. In addition, as a modified gravity theory, it needs 
to pass (3) the solar-system test, such as the constraints on the Brans-Dicke theory. The 
process of these three tests in the present paper are described in the following. 

(1) Cosmic Expansion History: For the background cosmic expansion there is a 
degeneracy between f(R) modified gravity and GR with dark energy, that is, one cannot 
distinguish between these two theories via the observations about the cosmic expansion 
history such as type la supernovae (SNIa) and baryon acoustic oscillations (BAO). This is 
because for any given expansion history H(t) one can always construct a function f(R) to 
generate the required expansion. The f(R) model so constructed is called "designer f(R)" 



22J. To fit the observational constraint on H(t), we take the designer f(R) approach to 



(re) construct f(R) for possible expansion histories, where f(R) is specified by H(t) and the 
initial condition, particularly the value of df jdR at some initial time, as denoted by f^. In 
the present paper we consider the expansion history H(t) given by the effective dark energy 
equation of state w e g as specified by the Chevallier-Polarski-Linder (CPL) parametrization 
23L l24j . wqpl = wo + w a z/(l + z) (where z is redshift), with the current observational 
constraints on w and w a [5]. In this case the designer f(R) is parametrized by three 
parameters, {wo,w a , fin}. We then test such designer f(R; {wo, w a , fm}) model via the 
cosmic-structure observations and the solar-system experiments. 

(2) Cosmic Structure Formation: For each designer f(R;{wo,w a ,fm}) obtained 
above we calculate the theoretical prediction of the cosmological perturbations, particularly 
and G e s/GN at late times and at subhorizon scales, where \I/ and $ are two metric 
potentials in the Newtonian gauge, and G e g and Gn are the effective gravitational coupling 
and the Newtonian constant, respectively. In GR these two ratios are both unity, while in 
f(R) modified gravity they generally change with time and scales. For the test we compare 
the theoretical prediction with the observational requirements from the cosmic structure 
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observations, thereby constraining f(R; {wo,w a , fm}), i.e., obtaining the constraints on the 
parameter space {wq, w a , fm}- In particular, we obtain the allowed range of fm for various 
{w ,w a } under consideration. 

(3) Solar-system Tests: The solar-system experiments provide stringent constraints 
on the Brans-Dicke theory with the chameleon mechanism, and accordingly tightly constrain 



the behavior of f(R) for R/Hq > 10 5 25l428j . To do the solar-system test, for each designer 
f(R; {wo, w a , fm}) we check whether it is consistent with this constraint on f(R/Hl > 10 5 ). 

Through the above three tests we find that only the models which closely mimic GR with a 
cosmological constant can survive these three tests. That is, among the f(R; {w , w a , fm}) 
models under consideration, only the case where f(R) — constant or {wo,w a , fm} very 
close to {—1,0,0} survives these three tests. For comparison, when we consider only the 
cosmological tests, i.e. the tests (1) and (2), for each {wo,w a } under consideration we find 
the allowed range of fm very narrow although existing. This presents a severe fine-tuning 
of the initial condition of f(R) for passing the cosmological tests. 

This paper is organized as follows. In Sec. HH we introduce the basics of the background 
evolution in the f(R) theory, the construction of the designer f(R) models, and the evolution 
of the cosmological perturbations with the late-time and the subhorizon approximations. In 
Sec. IIHI we introduce the observational constraints set by the observations about the cosmic 
expansion and the cosmic structures and by the solar-system experiments. In Sec. IIVI we 
show the results of the three tests. We then conclude in Sec. El 



II. BASICS OF f(R) MODIFIED GRAVITY 
A. Background evolution 

We start with the following modified Hilbert-Einstein action for the f(R) theory of mod- 
ified gravity. 1 

s = ^- 2 I d*xV=g [R + f(R)] + J2s a , (1) 

where f(R) is a general function of the Ricci scalar and represents the deviation from GR, 
and S a with a = r,m are the actions for radiation and matter (including dark matter and 

1 k 2 = 8irG 
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baryons), respectively. The variation of this action with respect to the metric tensor gives 
the gravitational field equations, 

GV + UrRulu - \f9^j + (9^D - V M V„) f R 

=«%-« 2 E T ^ ( 2 ) 

a 

where the d'Alembertian operator □ = V U V U , f R = df/dR, likewise = d 2 f/dR 2 , and 
the energy- momentum tensor of the cosmic fluid T$ = — (2/y/—g)(5S^/5g^ I/ ). 

Regarding the background evolution, we consider the flat Robertson- Walker metric for 
describing the background space-time, 

d s 2 = -dt 2 + a 2 {t)dx 2 , (3) 

and the perfect fluids for the energy contents, 

T fl L , = ^2diag(-p a ,p a ,p a ,p a ), (4) 

a 

where p a and p a are the energy density and pressure. We then obtain two evolution equations 
describing the cosmic expansion: 

H 2 + 1/ + (h 2 - V) f R + H f R = ^( Pm + Pr ), (5) 

- + \f - H 2 f R + ~Hfn + \f R = zf {pm + 2pr ) , (6) 
a o lib 

where an over-dot denotes the derivative with respect to the cosmic time t, the Hubble 
expansion rate H = a/a, the Ricci scalar R = 6H + 12H 2 , and we have used p m = and p r = 
pr/2>. Equation ([6]) shows that matter and radiation can simply contribute to deceleration, 
while the deviation from GR may induce accelerated expansion when its contribution therein, 
//6 — H 2 fji + Hfn/2 + fn/2, is negative and significant. In addition, from these equations 
one can see that the f(R) theory comes back to GR (with a cosmological constant A) at the 
background level when f(R) = constant = —2 A. 

The above two equations can be recast in the form in GR: 

K 2 

H 2 = y (p m + Pr + Peg) , (7) 
- = 7T (j>m + 2Pr + PeS + 3p e ff) • (8) 

a o 
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Here p e s and p e s are respectively the effective energy density and the effective pressure 
contributed from the deviation from GR, f(R), and are defined as 



Pes = 



PeS = 9 



lf-(H 2 + j)f R + 2Hf R + f, 



R 



(9) 
(10) 



This effective energy density satisfies the conservation equation, 



where 



PeS = -3H (1 + W e s) Peff • 



W cS = Peg/ PeS 



(11) 



(12) 



For a given expansion history H{t) and given energy densities, p m (t) and p r (t), Eq. 
([5]) gives a second-order differential equation for f(R) or f(t), with the understanding that 
fn — f JR an d fn = f / R — Rf / R 2 . Equivalently, for a given effective energy density p c g that 
together with p m and p r specifies the expansion history, Eq. ([9]) gives the same differential 
equation for f(R) 



f" 



H 2 



+ 



E R" 

H 1 

2E R' 



12E' 3E" 
+ 



JL 

H 2 o 

J_ 

6H 2 



cS 



0. 



E E 

where the prime denotes the derivative with respect to In a, and 
E(a) = H 2 /H 2 = n m a- 3 + fi r a -4 + E eS , 



(13) 



E eS (a) = pes/Pc = ^esexp <^ -3 



In a 



[1 + itfeff(a)] din a 



(14) 
(15) 



Here the density fraction Qi = pj(now) / p c (for i = m,r, eff), p c is the present critical density, 
and Q m + Q r + Q c g = 1 under our consideration of a flat universe. Accordingly, the function 
f(R) which satisfies this differential equation can generate the required expansion history. 
Since in general the solution of this differential equation exists (and is not unique), for 
any given expansion history there always exists a function f(R) which can generate that 
required expansion. As a result, we cannot distinguish between f(R) modified gravity and 
dark energy via the observations about the cosmic expansion. 
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B. Designer f(R) 



In the study of f(R) modified gravity, one may start with a carefully chosen function 
f(R)- On the other hand, one may instead reconstruct f(R) from a given expansion history 
(together with the information about the energy contents), i.e., considering the solutions / 
of Eq. (fl3|) for given cosmological parameters {H Q , Q m , Q r , fl e a} and given w e s(a). The f(R) 
so constructed is called "designer f(R)" 22]. This phenomenological approach makes sense 
particularly when one expects that the future cosmological observations will give precise and 
detailed information about the cosmic expansion. 

Solving Eq. ( II 3p requires initial conditions, with different setting of which one obtains dif- 
ferent solution f(R). Thus, the designer f(R) is in general a functional of {H , Q m , Q r , fl e s}, 
w e ff(a) and the initial conditions. In the present paper we will fix {Hq, Q m , Q r , fl e s}, while 
consider a variety of w^{a) and initial conditions for possible f(R). 

It has been shown 22] that the two initial conditions originally needed for solving the 
second-order differential equation ([13]) can be simplified to only one condition. This is 
done as follows. At the early times we require the deviation for GR be tiny in order to 
satisfy the stringent constraints about cosmic microwave background (CMB) and big-bang 
nucleosynthesis (BBN), and therefore p e ff "C p r ,m- With this approximation, one can obtain 
an analytic solution of Eq. (Tl3|) . which is the sum of a particular solution and a growing- mode 
and a decaying mode homogeneous solutions with two arbitrary constants to be determined 
by two initial conditions. Since the decaying mode entails a large deviation from GR at 
early times and cannot fit the above-mentioned requirement, one may exclude the decaying- 
mode homogeneous solution. Then only one arbitrary constant remains and only one initial 
condition is needed. The initial conditions for / and /' at an initial time a« can be specified 
as follows 



22|. 



-± A+a? + A p E eS ( ai ) , (16) 
f 

^2 ~ p+A+a^ - 3 [1 + WeffK)] A p E e $( ai ) , (17) 
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where 



—6c 

A, 



" ' -Sw'{a i )+9w 2 {a i ) + {lS-3b)w eS {a i ) + 9-3b + c , 



-b + Vb 2 -Ac 

p+ 



b 



2 

7 + 8r _ -3 
2(1 + r) ' ° = 2(1 + r) 

dgg 1 



The remaining one arbitrary constant is the amplitude of the growing-mode homogeneous 
solution, A + . Instead of specifying the value of A + , one may specify the more physical initial 
condition fn(ai) that gives /j and f[ through the relation fn = f'/R', Eqs. (|T6|) and (jTj 



C. Cosmological perturbations 

As pointed out at the end of Sec. Ill Al f(R) modified gravity and dark energy can lead 
to the same evolution of the background space-time, and therefore we cannot distinguish 
these two theories of cosmic acceleration via the observations about the cosmic expansion 
history. On the contrary, even for those leading to the same cosmic expansion, these two 
theories may have distinct effects on the cosmic structure formation, which accordingly gives 
an opportunity of breaking the degeneracy. 

For the cosmic structures at late times, the relevant physical quantities we consider here 
are the sub- horizon modes of the matter density perturbations 8 m , and the scalar metric 
perturbations in the conformal Newtonian gauge, \I> and $, defined by 

ds 2 = - [1 + 2*(x, t)} dt 2 + a 2 [1 - 2$(x, t)} dx 2 . (18) 

These two potentials respectively characterize the particle acceleration, 'x = — V\l/, and the 
gravitational potential, V 2 $ = AirG^p, in the Newtonian gravity. We consider the Fourier 
transformation of these three perturbed quantities, and in the remaining of the paper 5 m , 
\1/ and $ will denote the Fourier modes in fc-space with the understanding that they are 
functions of time t and wave-number k. 

With the sub-horizon and the late-time approximations, we have the following equations 



for the evolution of S m and the relations between 5 m , ^ and $ 



29) 



C+(2 + f)^-^^5 m ^0, (19) 



~ -4tt (a 2 /k 2 ) G cSPm 5 m , (20) 



(21) 



where the effective gravitational coupling strength G e g is given as 

(22) 



r> i 14- 4fc Jflfl 

<Jeff _ J- x ~ a 2 l+f R 



G N l + fRl + ^ R - 

Here Gn is the Newtonian gravitational constant measured in the solar system, which is 
very close to k 2 /87t because the deviation from GR in the solar system must be very tiny 
(if nonzero), i.e., <C 1. 

The two ratios, and G e s/GN, give the key to differentiating modified gravity from 
GR. In GR both of them are unity, while in f(R) modified gravity they are generally time- 
and scale-dependent. Equations fl2Tl) and ([22"]) give the theoretical prediction of these two 
ratios for a given f(R), which can be conveniently used for the comparison with the existing 
observational constraints to be presented in the next section. 



III. OBSERVATIONAL CONSTRAINTS 



In this section we present the observational constraints about (1) the cosmic expansion 
history, (2) the cosmic structure formation, and (3) the deviation from GR in the solar 
system. 



A. Cosmic expansion 

In regard of the effects on the cosmic expansion, a f(R) model is equivalent to a dark 
energy model with the effective equation of state w e g given by Eqs. ffXDT) and (|T2l) . 
Accordingly, the constraints on the equation of state of dark energy from the observations 
about the cosmic expansion also constrain f(R) or, more precisely, the designer f(R) models. 
In the present paper we consider the constraints [5( on a constant w e g and on the CPL 
parametrization i^cpl = Wq + w a z/(l + z): 

w eS = constant = -0.980 ± 0.053 (68% CL) . (23) 

Wq = -0.93 ± 0.13 , w a = -0.41+°;^ (68% CL) . (24) 
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We will consider the designer f(R) models which are constructed with respect to a variety 
of Wes following the above constraints. Such designer f(R) models by construction can pass 
the cosmic-expansion test. 



B. Cosmic structure formation 



Recently many works on the test of modified gravity via the observations about the cosmic 



structures have been done 



30h35] . In the present paper, for the test of f(R) modified gravity 



we will invoke the observational constraints on and G^/Gn obtained in 33] from the 
observational data about the integrated Sachs- Wolfe (ISW) effect, which is measured by 
cross-correlating the CMB with the galaxy data (the tracer of the large-scale structure) 
with the help of the type la supernova (SN la) data for the information about the cosmic 
expansion. 



In [33] the cosmic background expansion is fitted with the ACDM model, and the cosmic 
perturbations and G^/Gn are parametrized by some fitting formulae with only one 
free parameter Ai that corresponds to the mass scale of the scalaron. 2 The constraint on Ai 



obtained in |33|| is: 

Ai < 1900Mpc//i at 95% c.l. , (25) 

which then gives the allowed ranges of and G e s/G^ for various time and scales k. In the 
present paper, for the test of f(R) modified gravity we will simply consider the constraint 
for the present time (a = 1) and the scale k = 0.01/iMpc -1 : 

1< */$ < 1.996, (26) 
1 < G cG /G N < 1.403. (27) 

This constraint is consistent with GR where both ratios are unity. We note that according 
to our experience the test of f(R) modified gravity is not sensitive to the choice of the scales. 

We expect the future observations will put more stringent bounds to and G e ^/G^. 
For instance, Guzik et al. {32] use lensing power spectra, galaxy-lensing and galaxy-velocity 
cross spectra to estimate the ability of the future surveys to test gravity theories. With 



2 The scalaron is a scalar degree of freedom introduced when the f(R) theory is transformed to the scalar- 
tensor theory. 
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ACDM as the fiducial model they give the following narrower allowed range: 

0.83 < */$ < 1.25 , 
0.84 < G cS /G N < 1.18. 



(28) 
(29) 



C. Solar-system constraints 



The solar-system experiments severely constrain modified gravity and the deviation from 
GR, including the f{R) theory, in the regime with the Ricci scalar R/Hq > 10 5 . From 



these local experiments Gu and Lin 
the chameleon mechanism: 



28| have obtained a general constraint on f(R) with 



< RfRR < 0.4 as R/Hl > 10 5 
10- 16 < f R <0 as R/Hl ~10 £ 



(30) 
(31) 



IV. THREE TESTS OF f(R) MODIFIED GRAVITY 

In this section we test the f(R) theory of modified gravity via the observations about 
the physics at three different scales, ranging from the largest scale to the local scale: (1) 
the cosmic background expansion history, (2) the cosmic structure formation, and (3) the 
solar-system test of gravity. 



A. Cosmic-expansion test 

To fit the requirement from the observations about the cosmic expansion history, we 
consider the designer f(R) constructed with respect to the given cosmological parameters 
{fl m , Q r , ^eff} consistent with observations, the effective equation of state w e s satisfying the 
constraint in Eqs. fl23l) or (|24"|) . and various initial conditions fm = /j?(aj). In particular, we 
set the parameters Q m = 0.27, Q r = 8.4 x 10 -5 , f2 e fr = 1 — Q m — &>r for a flat universe and 
the initial time = 10~ 8 , and consider a variety of constant w e g and the CPL-parametrized 
effective equations of state, w e g(z) = w + w a z/ (1 + z), as follows. 

• For the const ant-u^g- case, following the constraint in Eq. (|2"3"|) we set w a = and 
choose 6 different values for u> : 
w = -1.03, -1.00, -0.99, -0.97, -0.95, -0.93. 
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• For the CPL-w e s case, following the constraint in Eq. f )24|) we choose four wq and eight 
w = -1.04, -1.00, -0.93, -0.85, 

w a = -1.10, -0.72, -0.41, -0.22, 0.00, 0.15, 0.25, 0.30, 
which together give 32 different u> e ff- 

As an example for demonstration, Fig. [I] shows the designer f(R) with respect to w c s = — 1: 
five f(R) corresponding to five initial conditions fm- 

We then test these designer models, f(R; {w , w a , fm}), via the cosmic-structure obser- 
vations and the solar-system experiments, as to be shown in the next two subsections. 

B. Cosmic-structure test 

To test the designer f(R) models via the observations about the cosmic structure forma- 
tion, for each model we check whether the theoretical prediction of the two ratios, and 
Ges/GN, m Eqs. ( 12T]) and ( 122]) can fit the observational constraint in Eqs. ( )26|) and ( )27|) for 
the present time (a = 1) and the scale k = 0.01/iMpc -1 , which is the only case (time and 
scale) we consider in the remaining of the paper. We note that the test of f(R) modified 
gravity is not sensitive to the choice of the scales. 

n 

This observational constraint is deduced from [33| where the cosmic expansion history 
is fitted with the ACDM model. Accordingly, strictly speaking, it can only constrain the 
designer f(R) with respect to w e ff = —1. Nevertheless, here we also apply this constraint to 
other models in order to manifest how various designer f(R) models might be constrained 
by the current cosmological observations. 

1. 1 < /$ < 1.996 for {a = 1, k = O.Ol/iMpfT 1 } 

For each w e ff under consideration we explore numerically the relation of the ratio ^//^(a = 
1, k = O.Ol/iMpc -1 ) to the initial condition fm. As an example for demonstration, in Fig. [2] 
we show the relation of to fm for the designer f(R) model with respect to w e s = — 1, 
as compared to the observational constraint denoted by the grey region. In addition, we 
show in Fig. [3] the relation for various constant w e ff, and in Fig. H] that for the best fit of the 
CPL-parametrized w e ff- 
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We find this relation similar for different w e s: = 2 for almost all fm, while it changes 
violently and goes to ±00 within a narrow range of fm, only in which the designer models 
may fit the observational constraint in Eq. ( I26p . As a result, one needs to fine-tune the 
initial condition in order to have a designer f(R) model consistent with the observational 
results about the cosmic structure formation. 

We illustrate the narrow viable range of the initial conditions fm in Figs. [5] and |6] for 
the designer f(R) models with respect to the constant w e s and the CPL-parametrized w e g, 
respectively. We write the viable range as fm = fm ± Afm, and present in these plots the 
central value ±fm (denoted by + and x) and the width 2Af Ri (denoted by *) of the viable 
range. 

As have been pointed out above, most of the designer f(R) models predict = 2, 

which is around the margin of the current constraint fl26l) . We expect the future cosmological 
observations will give more detailed and precise information about cosmic structures, and 
therefore a sharper bound on If the future upper bound of \&/$ is significantly smaller 

than 2, the f(R) modified gravity will be even more fine-tuned and therefore disfavored. On 
the other hand, if the future observations give a bound excluding unity (the GR prediction), 
GR will be ruled out. 

2. 1 < G eff /G N < 1.403 for {a = 1, k = O.OlhMpc- 1 } 

For each w e g under consideration we explore numerically the relation of the ratio 
G e ff/G ! 7v(a = 1, k = 0.01/iMpc^ 1 ) to the initial condition fj$. For demonstration, we show 
in Fig. [7] the relation for the designer f(R) model with respect to w c g = — 1, and in Fig. [8] 
that for the best fit of the CPL-parametrized w e g, as compared to the observational con- 
straint denoted by the grey region. We find this relation similar for different w e Q\ G e s_/G^ 
smoothly changes with fm for most of the initial conditions, but varies violently and goes 
to ±00 within a narrow range of fm, around which the designer models can fit the observa- 
tional constraint in Eq. (I27j) . As shown in these two figures, there are two separated viable 
regions of the initial conditions, one narrow region in the vicinity of the singularity and the 
other wider region extended to the initial conditions away from the singularity. As a result, 
the test with the ratio G e g/GN gives a weaker constraint on the designer f(R) models under 
consideration. 
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C. Solar-system test 

To perform the solar-system test of the designer f(R) models, for each of the models 
under consideration we check whether the requirement in Eqs. (13"U1) and (T5TT) is satisfied. 
We examine the designer models constructed with respect to the w e s listed in Sec. IIV Al 
and a wide range of the initial conditions fm. As a result, we find the constraint from the 
solar-system test much more stringent: only the models with w e ff = —1 and fm = 0, i.e. 
f(R) = constant, can fit the constraint. These viable models are very close to the case of 
GR with a cosmo logical constant, and therefore cannot be differentiated from the ACDM 
model, both cosmologically and locally. 

To manifest how close to GR the viable models need to be, we survey in details the small 
portion of the parameter space (w e g, fm) around the GR point (—1, 0) for constant w e s. We 
present the viable region in Fig. [9j This figure shows that even the cases barely different 
from GR are not viable. Thus, among the designer f(R) models under consideration, only 
the models closely mimicking GR can pass the solar-system test. Furthermore, we note that 
the viable models are so close to GR that in these models not only the cosmic expansion 
is nearly ACDM, but also the key quantities for distinguishing GR from modified gravity 
in cosmic structures, and G c s/Gn, agree with the GR prediction (unity) better than 
10 -6 . Accordingly, these models cannot be differentiated from GR via the cosmological 
observations, including the observations about the cosmic structure formation and those 
about the cosmic expansion history. As a result, the solar-system test of gravity rules out 
the frequently studied w e s = — 1 designer f(R) models which are cosmologically distinct 
from ACDM. 

V. CONCLUSION AND DISCUSSION 

In this paper we have investigated three tests of the f(R) theory of modified gravity, 
including two cosmological tests and one local test. We have invoked the constraints from 
the observations of (1) the cosmic expansion history and (2) the cosmic structure formation 
and from (3) the solar-system experiments. We found the constraint from the solar-system 
test particularly stringent. 

To pass the cosmic-expansion test we have considered the designer f(R) models con- 
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structed with respect to the constrained effective CPL equation of state w e s(z) = wq + 
w a z/(l + z) and various initial conditions fm = /#(a = 10~ 8 ), and checked whether they 
pass the other two tests. We have presented the constraint on from these two tests for 
each w e g(z). We conclude that the solar-system test rules out all of the designer f(R) models 
constructed with respect to wcpl, except those closely mimicking GR with a cosmological 
constant, i.e., the ACDM model. In addition, even passing the cosmic-structure test alone 
would require a fine-tuning of the initial condition of f(R)- 

We have particularly explored the viable region in the designer-/ (R) parameter space 
(w e ff, fm) around the GR point (—1,0) for constant w e g, and found it extremely small. Even 
a tiny deviation from the GR point is ruled out by the solar-system test. The designer models 
within this viable region are cosmologically indistinguishable from the ACDM model. More 
precisely, their prediction of the cosmic expansion and the cosmic structure formation agrees 
with ACDM better than 10~ 6 . As a result, the frequently studied w e s = —1 designer f(R) 
models that are cosmologically different from ACDM have been ruled out by the solar-system 
test. 

We note that despite the stringent constraint from the solar-system test, it is still possible 
to construct a viable f(R) model that is observationally different from GR. As pointed out 
by Gu and Lin in [28j, the solar-system test gives a severe constraint [Eqs. (130]) and (T3~TI) ] on 



the behavior of the function f(R) for R/Hq around or larger than 10 5 , which corresponds 
to the time around or before z ~ 50. Thus, roughly speaking, the solar-system test and the 
observations about CMB and BBN require the deviation from GR be extremely small, i.e. 
f(R) — constant, from the early time to z ~ 50, while a significant deviation from GR is 
allowed in the present epoch. Accordingly, a viable designer f(R) model can be constructed 
with respect to the w e g which is close to —1 in the past but significantly differs from — 1 
recently. These designer f(R) models are worth consideration in the quest for a modified 
gravity model that is cosmologically distinct from GR with a cosmological constant. 
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FIG. 1: Five different designer f(R) w.r.t. w e s = —1 for different initial conditions fj^. 
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FIG. 2: The relation of to fm for the w e g = —1 designer f(R) model, as presented by the 
solid curve. The grey region shows the observational constraint of and the dot denotes the 

GR prediction, \I>/<3> = 1. 



18 



© 1.8 





observational constraint 

"Iff— i -oo 

w eff = -1.03 

w efr --0.99 

w eff = -0.97 

w eff = -0.95 

w cll = -0.93 , 
















■ 





















-18 -16 -14 -12 -10 -8 -6 -4 -2 



FIG. 3: The relation of to A + for the designer f(R) models w.r.t. u> eff = -0.93, -0.95, -0.97, 
—0.99, —1.03, —1.00, as presented by the solid curves from the left to the right, respectively. The 
constant A + is related to the the initial condition of f{R) through Eqs. (|16|) and (|17|). The grey 
region shows the observational constraint of Vl//<&. 
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FIG. 4: The relation of Vl//$ to fm for the designer f(R) models w.r.t. the best fit of the CPL- 
parametrized u) e fj, i.e., wo = —0.93 and w a = —0.41. The grey region shows the observational 
constraint of ^/$. 
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FIG. 6: The viable initial conditions fja = fni±Afjn for the CPL-^eff designer f(R) models. The 
central value ±fm is denoted by + or x, and the width 2Afja of the viable range denoted by *. 
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FIG. 7: The relation of G^q/Gn to fm for the w e g = —1 designer f(R) model, as presented by the 
solid curve. The grey region shows the observational constraint of G e s/GN, and the dot denotes 
the GR prediction, G e s/G]y = 1. 




FIG. 8: The relation of G e s/G^ to fm for the designer f(R) models w.r.t. the best fit of the 
CPL-parametrized w e s, i.e., wo = —0.93 and w a = —0.41. The grey region shows the observational 
constraint of G^/Gn- 
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FIG. 9: The viable f(R) models in the parameter space (w e g, fm) for constant w e g under the solar- 
system constraint. The viable region is presented by the red contours in the four plots together. 
The GR point (—1, 0) is between the four plots. 
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